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The iconic Schrodinger's cat state describes a system that may be in a superposition of two macro- 
scopically distinct states, for example two clearly separated oscillator coherent states. Quite apart 
from their role in understanding the quantum classical boundary, such states have been suggested 
as offering a quantum advantage for quantum metrology, quantum communication and quantum 
computation. As is well known these applications have to face the difficulty that the irreversible 
interaction with an environment causes the superposition to rapidly evolve to a mixture of the com- 
ponent states in the case that the environment is not monitored. Here we show that by engineering 
the interaction with the environment there exists a large class of systems that can evolve irreversibly 
to a cat state. To be precise we show that it is possible to engineer an irreversible process so that 
the steady state is close to a pure Schrodinger's cat state by using double well systems and an envi- 
ronment comprising two-photon (or phonon) absorbers. We also show that it should be possible to 
prolong the lifetime of a Schrodinger's cat state exposed to the destructive effects of a conventional 
single-photon decohering environment. Our protocol should make it easier to prepare and maintain 
Schrodinger cat states which would be useful in applications of quantum metrology and information 
processing as well as being of interest to those probing the quantum to classical transition. 



The development of many quantum technologies de- 
pends on an ability to engineer strongly non classical 
states. Such states take the form of either highly entan- 
gled states of distinct degrees of freedom or a quantum 
coherent superposition of macroscopically distinct states 
in a single degree of freedom 1 , known as Schrodinger's 
cat states (after a well known thought experiment 2 ). It is 
these cat states that we consider in this letter. There has 
been great progress in the production of such states as 
well as experimentally reconstructing such states through 
a series of measurements in a process known of as quan- 
tum state tomography 3 8 . These developments are of 
great importance as, in addition to their curious nature, 
Schrodinger cat states can be used as a resource for devel- 
oping technologies such as quantum computing 9 ' 10 , quan- 
tum communication 11 ' 12 and quantum metrology 13 15 . 
The main obstacle to deploying cat states in such appli- 
cations is their fragility as they are destroyed by noise in 
a process termed environmental decoherence. A careful 
consideration of optical cat states shows that this deco- 
herence may be interpreted as due to Poisson distributed 
jumps between even and odd cat states whenever a single 
photon is lost 16 18 . Their production and maintenance 
requires very precise quantum control as well as low dis- 
sipation. In this work we propose a protocol for double 
well systems to create Schrodinger cat states that actu- 
ally uses the non-controllable, non-unitary interaction of 
the system with a special kind of environment to create 
Schrodinger's cat states. To be specific, we have found 
that for a simple double-well system system interacting 
with an environment comprising a bath of two-photon 
absorbers, for certain initial states, the system relaxed 
to a steady state which is close to a pure Schrodinger cat 



state. Such an environment when paired with a paramet- 
ric photon pump is known to exhibit many interesting 
effects in quantum optical systems, from cats to quan- 
tum statistics 19 ' 20 . Two-photon absorption has also be 
suggested as a powerful resource for quantum comput- 
ing application 21 . Two-photon decay preserves parity 
and enables the system to relax to a steady state with 
same parity as the initial state. Our model is simpler 
than and different from other driven dissipative bistable 
systems (for example, the coherently driven optical cav- 
ity containing a Kerr medium 22 , the driven Duffing me- 
chanical resonator 23 , tapered optical fibers 24 and photon 
pumps 19 ' 20 ), as we do not include driving on either the 
cavity resonance or the coordinate degree of freedom. 
Our proposal opens up new opportunities for explor- 
ing quantum phenomena from the micro to macroscopic 
level and in fields as diverse as quantum optics 25 , Boes- 
Einstein condensates 26 , quantum electronics 27 and nano- 
mechanics 28 (for which multi-phonon relaxation has al- 
ready been proposed 29 ) or any other system in which it 
is possible to generate a double well potential. 

For the results presented in this paper we have used 
as an example system a superconducting quantum in- 
terference device (SQUID) ring. Our reason for choos- 
ing SQUIDs is that these devices are routinely fabricated 
and their theory is very well understood. We note that 
we have investigated a number of other systems (but do 
not include results here) and our analysis indicates that 
the key feature of the ring is that it can be made to form 
a double well potential. Moreover, nonlinear systems de- 
rived from the Josepheson junction in circuit QED ex- 
hibit multi photon resonance when driven by an external 
field 30 and thus we expect two-photon decay to be present 
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in such systems. The real difficulty is making it dominate 
over single photon effects. We will return to this later. 
Beyond these considerations we believe there is nothing 
particularly special about the exact form of the potential 
needed to realise our protocol. Subject to being able to 
engineer an appropriate dissapative channel we therefore 
believe that the methodology that we propose for gener- 
ating cat states will, as previously mentioned, find wide 
application. The potential energy of the SQUID compris- 
ing a thick superconducting ring enclosing a Josephson 
junction weak link takes the form of a harmonic oscillator 
perturbed by a cosine 

U($ x ) -cos 2tt— - 
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where the coordinate <1> is the total magnetic flux in the 
ring and <I>o — h/2eis the superconducting flux quantum. 
We have chosen example circuit parameters that are in- 
line with modern fabrication techniques and suited to 
experimental realisations: A = 3 x 10 _10 H for the ring's 
inductance and I c = 2pA as the critical current of the 
weak link (although not in the above formula we also 
chose a capacitance C — 5 x 10 _15 F). We set the ex- 
ternally applied magnetic flux $ x = 0.5$o so that the 
ring's potential forms a degenerate double well. It is 
also convenient to introduce the bosonic annihilation a, 

and creation operators <I> = J ' 2 cu LC ( a wnere 

cjlc — 1/VAC. In Fig. 1 we show the potential energy 
of the ring as well as the energy of the ring's station- 
ary states. It is worth noting that the ground state and 
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FIG. 1: Stationary state energy levels: The potential 
energy of the ring (black) as well as the energy of the rings 
stationary states (blue). Parameters used here and through- 
out the paper are inductance A = 3 x 10 _10 H, capacitance 
C = 5 x 10~ 15 H, critical current of the weak link I c = 2pA 
and externally applied magnetic flux <fr x = 0.5$o- Note that 
we have exaggerated the energy difference between the ground 
and first excited states as well as stationary states two and 
three in oder to make the different energies visible on this 
plot. 



first excited state approximate, respectively, symmetric 
and anti-symmetric superpositions of two coherent states 
centred at the bottom of each well. These two states have 
very nearly the same energy and the difference in their 
energy has been exaggerated in this plot (as have those 
for the second and third excited states). 

We model the effect of the environment on the system 
using the master equation in the Lindblad form 31 
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where p is the density matrix describing the state of the 
system (initially p — \ijj(t = 0)) (rf(t = 0)|) and H is the 
system's Hamiltonian. The non-unitary effect of the en- 
vironment on the system is contained in the Lindbald op- 
erators Lj with each describing a possible environment. 
For example the usual Ohmic (i.e. analogous to fric- 
tion proportional to velocity) bath, at zero temperature, 
would be described by a Lindblad proportional to the an- 
nihilation operator. For an undriven system the master 
equation has steady state solution that, in the presence of 
an environment, is usually a density operator in a mixed 
state. In certain circumstances, at zero temperature, 
these solutions may be pure states such as the vacuum 
state of the harmonic oscillator. In these circumstances 
the solutions will not exhibit features such as superposi- 
tions of macroscopically distinct states and are relatively 
uninteresting. It is precisely this process where the envi- 
ronment essentially removes the system's quantum coher- 
ence from de-localised, or more generally non-Gaussian, 
states that is known of as environmental decoherence. 
The density matrix for a decohered system without these 
quantum correlations represents a statistical mixture of 
possible states of the system and, for a single quantum 
object, can be directly compared with classical probabil- 
ity density distributions 32 . It should be noted however 
that there are driven dissipative systems, for example dis- 
persive bistability, for which the steady state is a mixed 
state with a considerable amount of quantum coherence 
in the limit of large Kerr nonlinearity 16 ' 17 ' 33 . 

We found very different behaviour if one chooses a 
different environment comprising two-photon absorbers, 
described by a Lindblad proportional to the square of 
the annihilation operator. In Fig. 2 we show the en- 
ergy expectation values and von-Neumann entropy, S = 
— Tr[p In p] as functions of time for solutions of the master 
equation for the ring in the presence of such an environ- 
ment. We used as initial conditions the first twenty en- 
ergy eigenstates of the ring Hamiltonian. In these plots 
the energy behaves just as one would expect the energy 
of an undriven open quantum system to do - it settles to 
a single value. When one inspects the dynamics of the 
entropy however the story is quite different. One usually 
expects the entropy to grow from zero to some asymptotic 
value as the system evolves into a mixed state. While we 
see that this is the initial behaviour the entropy does not 
monotonically increase, instead it decreases until the en- 
tropy is nearly negligible. It appears that the the system 
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has to a significant extent recohered and the final density 
matrix is very nearly that of a pure state. While this is 
not the usual behaviour of an open quantum system it 
is in-line with our expectations of an environment that 
"decoheres" a system to an almost pure state that is a 
very good approximation to a Schrodinger cat state 19 . 

In order to demonstrate that the system does indeed 
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FIG. 2: Effect of decoherence on energy and entropy 

We show the dynamical evolution of the ring's energy and 
entropy using each of its first twenty stationary states as ini- 
tial conditions. The dynamics have been found by solving 
the master equation for the ring in the presence of a bath of 
two-photon absorbers (with L = y/0.2a 2 ). We have provided 
insets for increased resolution of the system's initial dynamics. 
The top plot shows the dynamics of the ring's total energy. 
As expected for an open quantum system of this kind the ring 
can be seen to decohere to one energy, a little above that of 
the ground state. The bottom plot shows the dynamics of the 
von-Neuman entropy for the ring. In each case the initial en- 
tropy is zero as the system starts in a pure state. The entropy 
grows before dropping off to a low value indicating that the 
systems steady state solution is very nearly a pure state. 



decay to a Schrodinger cat state we will make use of 
the Wigner function. These pseudo probability density 
functions in phase space have been of great utility in 
demonstrating that some quantum states are Schrodinger 
cats 3 . The Wi emer function is 

W(*, Q) = ^J(* + CH* " exp (-^p) dC 

where Q is the charge variable that is conjugate to the 
magnetic flux <1>. In Fig. 3 we show three Wigner func- 
tions. Fig. 3a shows the initial state and is a coherent 
state centred at the origin. This is clearly recognisable 
as the expected Gaussian bell shape associated with co- 
herent states. We have solved the master equation for 
the ring in a lossy bath, with a Lindblad of L = \/(L2a 
and allowed the system to reach its steady state to obtain 
Fig. 3b. This is the Wigner function of a statistical mix- 
ture of two macroscopically distinct states and is in-line 
with expectations of the effect of a decohering environ- 
ment on such a device 34 . In Fig. 3c we show the Wigner 
function that we obtain by solving the master equation, 
as for (b), but replacing the damping term with a bath of 
two-photon absorbers, with L = . We notice two 

things: firstly that the state has rotated which we believe 
to be a consequence of a squeezing action associated with 
the bath and secondly that there are interference terms 
between the distinct states of the system. These interfer- 
ence terms, indicating quantum coherence, confirm that 
this state state is indeed a very good approximation to 
a Schrodinger cat. In oder to examine quantitatively the 
emergence of this cat from the initial coherent state we 
introduce, following 35 ' 36 , a measure of how de-localised 
the system is in phase space that is the integral of nega- 
tive parts of the Wigner function 



N(p) 
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In absolute terms this is a useful measure, but when we 
know (by inspecting the Wigner function) that the states 
we are examining are cat-like a more useful measure may 
well be a relative cattiness to some reference Schrodinger 
cat state. Hence we define: 



Cat(p,p re f) 



N(Pref) 



which quantifies the ratio of the de-localisation of one cat 
state against a reference cat and enables us to quantify 
if one is more [Cat(p, p re f) > 1], less [Cat(p, p re f) < 1] or 
just as [Cat(p, p re f) = 1] catty than the other. In Fig. 4 
we show the dynamics of this quantity for comparison 
with the results presented in Fig. 3 using as a reference 
state p re f the final cat state shown in Fig. 3c. Here we 
can clearly see that the cattiness of the system subject to 
an environment of two-photon absorbers monotonically 
increases and asymptotically converges to a steady state. 

It is interesting to consider what would happen to a 
ring that was initially in a Schrodinger cat state under 
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FIG. 3: Cooling for a cat In this figure we show, by making use of Wigner functions, the effect of two different environments 
on a ring prepared in a coherent state biased at zero flux. Each graph contains a top down view with a three dimensional plot 
of the function as a not to fixed scale inset. The graphs show a, the initial state which takes the form of a Gaussian bell, b, 
the steady state solution to the master equation under the influence of a conventional decohering environment comprising a 
lossy bath (with a Lindblad proportional to the annihilation operator L = y/0/2a). The ring has decohered to two distinct 
macroscopic states we do not see the interference terms between them that are characteristic of a Schrodinger cat state. We have 
instead a statistical mixture, the usual and expected result 34 , c, the steady state solution to the master equation for the ring 
coupled to a bath of two-photon absorbers (with a Lindblad proportional to the square of annihilation operator L = \/0^2a 2 ). 
In this case the ring has decohered to a superposition of two macroscopic states but now there are interference terms between 
these states indicating quantum coherence - the signature of a Schrodinger cat state. 



the influence of a bath of two-photon absorbers. For sys- 
tems with deep enough double well potentials such as the 
one considered here the ground and first excited energy 
eigenstates are both Schrodinger cats. The ground state 
is, to good approximation, an even superposition of two 
macroscopically distinct coherent states while the first 
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FIG. 4: Relative cattiness We show the cattiness measure 
Cat(p, p re f) for the the dynamics leading to Fig. 3b in red 
and to Fig. 3c in green. Here we have used as a reference 
state p re f the final cat state shown in Fig. 3c. For reference 
later we have also included the dynamics of Cat(p, p re f) for 
an environment of two-photon absorbers and damping. 



excited state is an odd superposition as can be seen from 
their Wigner functions in Fig. 5a and c respectively. The 
even and odd nature of these superpositions is reflected 
in the Wigner function by the phase of the interference 
terms between the two Gaussian's of the cat. It is known 
that such states would decohere under the environment 
of a lossy bath to a statistical mixture 34 . The dynamics 
of the system coupled to an environment comprising a 
bath of two-photon absorbers are, once more, found by 
solving the master equation with an L = v / 0^2« 2 , until an 
approximate steady state is reached. The Wigner func- 
tion of these states is then shown with Fig. 5a evolving 
to b and Fig. 5c to d. We observe that the phase in the 
final cat reflects that of the initial cat and the system 
has not simply decohered to the same steady state. The 
environment thus seems to preserve some of the symme- 
try of the initial state. We have checked the first twenty 
stationary states all of which decay to one of these cats 
or the other. Moreover, the pattern that was observed 
from the ground and first excited state persists and all 
even and odd states seem to evolve to cats of the same 
form as those shown Fig. 5b and Fig. 5d that are out of 
phase with each other. 

Our protocol seems all very well and good but an en- 
vironment of two-photon absorbers is very special. It 
would be hard to construct such an environment with- 
out having any other source of decoherence present. We 
therefore need to verify that the effects of a two-photon 
absorbing environment cannot be completely destroyed 
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FIG. 5: Preserving a cat Here we look at the ring initially 
in either its ground or first excited stationary state. As can be 
seen from their Wigner functions, plots a, and c, respectively, 
these take the form of Schrodinger cat states. The ground 
state is, to good approximation, an even superposition of two 
macroscopically distinct coherent states while the first excited 
state is an odd superposition. In terms of the Wigner func- 
tions this is reflected in the phase of the interference terms 
between the two Gauusian's of the cat. The effect of evolving 
the system in the presence of a bath of two-photon absorbers 
(L = \^0.2a 2 ) is then shown with a evolving to b and c to d. 
We observe that the phase in the final cat reflects that of the 
initial cat and the system has not simply decohered to the 
same steady state. 



by the presence of a more traditional environment such 
as a lossy bath. In Fig. 6 we show the results of just such 
a check. For each plot the system's initial state was the 
ring's ground energy eigenstate as shown in Fig. 5a. In 
Fig. 6a we show the effect of a lossy bath. We solve the 
master equation with a Lindblad L = v / 0.02a and allow 
the system to evolve until it has just decohered to a sta- 
tistical mixture and we have plotted the Wigner function 
at this point in time. We use this run as a benchmark 
for computing the next two cases which show the Wigner 
function solutions of the master equation computed over 
the same interval. In Fig. 6b we show the effect of a 
two-photon absorbing bath once more "decohering" to 
a Schrodinger cat state [L — ^/0^2a 2 ). And in Fig. 6c 
we apply both the lossy bath of in Fig. 6a and a two- 
photon absorbing environment of in Fig. 6b to the ring 
(Li = and L 2 = ^/^2a 2 ). We see that in this fig- 

ure there remain residual Schrodinger cat state features 
in the Wigner function. Hence, it seems that not only 
does a bath of two-photon absorbers create Schrodinger 
cat states, it also enables Schrodinger cat states to be 
more resilient to other forms of decoherence. In other 



words the presence of an environment of two-photon ab- 
sorbers seems to be prolonging the life of a damped cat. 
In Fig. 6d we quantify the cattiness using Cat(p, p re f) 
using the initial stationary state as shown in Fig. 5a as 
the reference cat. For the three environments considered 
here we find that for the system's later evolution the en- 
vironment of two-photon absorbers does indeed prolong 
the lifetime of the initial cat even in the presence of a 
lossy bath. We note that we obtain an almost identi- 
cal set of results if we start the system off in a coherent 
state centred at the origin (as in Fig. 3a). We chose to 
use the ring's ground state as, in our view, we obtained 
a more instructive plot of the states cattiness from the 
systems dynamics. For a direct comparison of the dy- 
namics of Cat(p, p re f) for these two initial conditions we 
now note that the dashed line shown in Fig. 4 was found 
for a lossy bath and a two-photon absorbing environment 
with L\ — \/0.02a and L 2 = V(L2a 2 . The green and blue 
lines of Fig. 4 and Fig. 6d are directly comparable. The 
idea that the presence of a two-photon absorbing envi- 
ronment can be used to extend the lifetime (and also 
generate) Schrodinger cat states holds equally well for 
two very different initial conditions. 

In order to make our above discussion a reality we 
need to engineer a dissipative quantum channel that 
acts as a two-photon absorber. Here we suggest a con- 
crete realisation that, whilst not perfect, still retains 
the key feature of environmentally induced "decoher- 
ence" to a Schrodinger cat state. Our proposal makes 
use of non-linearly coupled electromagnetic fields and 
SQUIDs. Such quantum electrodynamic circuits have 
already been investigated in the context of weak non- 
demolition measurement 37 ' 38 . One example comprises 
two microwave superconducting resonators coupled via 
a SQUID which in addition to a cross Kerr effect also 
manifests two photon conversion terms if the cavities are 
resonant 37 . Such systems can be quantised 39 42 and with 
a suitable arrangement and choice of circuit parameters 
can be reduced 43 ,43,44 to the form of a double well sys- 
tem subject to a two-photon absorbing environment (see 
supplementary material for details). Unavoidably, this 
process also brings with it an additional dephasing term, 
that adds to the master equation another Lindblad pro- 
portional to a^a. Nevertheless, we can report that whilst 
the dephasing term smears out the Gaussian peaks in 
the cat the interference terms in the Wigner function 
representing quantum coherence between the cat states 
remains strong. The fact that this dephasing term pre- 
serves parity is once more the key factor in ensuring the 
steady state of our engineered dissipative channel is still 
a Schrodinger cat state. Our proposal could lead to an 
initial realisation of a two-photon absorbing environment 
and concomitant interesting effects. The engineering of 
improved dissipative channels, without additional and 
unwanted decoherence effects, remains an open and in- 
teresting problem. 

There are two phenomena that embody quantum me- 
chanics, namely entanglement and the Schrodinger's cat 
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FIG. 6: A stubborn cat: combatting the effect of other forms of decoherence For each plot the system was initialised 
in its ground state of the ring as in Fig. 5a. In these plots we show a, the effect of a lossy bath on the state reducing a typical 
plot of a cat that has just decohered to a statistical mixture - setting the time that we use to sample the other two plots 
of this figure. (L = y / 0.02a). b, the effect of a two-photon absorbing bath showing decoherence to a Schrodinger cat state 
(L = V0.2a 2 ) c, the effect of both a lossy bath and a two-photon absorbing bath on the state. Notice that there are still 
signatures of a cat state unlike for the lossy bath alone - the environment of two-photon absorbers seems to be prolonging the 
life of the cat (L\ = V0.02a and L2 = V0^2a 2 ) and d, we show the cattiness Cat(p, p re f) for these three environments as a 
function of time (we have used the initial stationary state as shown in Fig. 5a as the reference cat in this case). We see for the 
system's later evolution the environment of two-photon absorbers does indeed prolong the lifetime of the initial cat even in the 
presence of a lossy bath. 



thought experiment 2 . The latter was proposed to high- 
light the difficulties we have connecting quantum me- 
chanics with everyday experience it neatly demonstrates 
the problems of understanding the emergence of the clas- 
sical world from quantum theory and the measurement of 
quantum systems. Schrodinger 's cat has become the icon 
of the subject and evolved to have a well defined meaning. 
It is an accepted explanation within the popular litera- 
ture that the reason the original thought experiment does 
not translate into reality (if conducted with a real cat in a 
box etc.) is that the environment to the radiation source 
(which included the cat itself) deletes the quantumness 
connecting the two states in a process known as decore- 
hence. As such environmental decoherence is something 
that many deem to be a crucial element in the the quan- 
tum to classical transition 32 ' 45 48 . We have shown that 
some environments may have a dramatically different 
effect on double well systems producing very quantum 
states as a result of "decoherence" . It may well be that 



system and environment such as the one we have used 
here could play an interesting role in quantum mechan- 
ically enhanced metrology probing foundational aspects 
of quantum mechanics associated with realising macro- 
scopic quantum phenomena and the quantum to classical 
transition. Furthermore, it is known that open quantum 
systems can be used to model the measurement process. 
Although it is beyond the scope of the current paper, we 
conjecture that it may well be possible to make use of an 
environment to measure a system into a Schrodinger cat 
state. 
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We show that a microwave superconducting cavity can be engineered to have a dominant two 
photon decay term using two cavities coupled by a SQUID. 

PACS numbers: 



There are a number of models [1, 2] whereby two mi- 
crowave superconducting cavities can be nonlinerally 
coupled using SQUIDs. We will base our discussion on 
Kumar and Divincenzo[l]. In that model, the hamilto- 
nian describing two microwave cavities, a probe (p) cavity 
and a signal (s) cavity, coupled with a SQUID is 

H = E Cp nl + E Lp 4> 2 v + E Cs n 2 s + E Ls 4> 2 s 
+ A [Ei p $ cos 4 P + Ei 8 <ft sin 4 /3 
+ 6El p Eljy s cos 2 p sin 2 p] (1) 

where n a ,(j) a are the standard charge and phase conju- 
gate variables describing the collective electrical degree 
of freedom in each cavity and A = 167t 2 Li/$q with 
L\ defined as the coefficient of the leading non-linear 
current term of the SQUID inductance. We will set 
cos 2 P = sin 2 P = 1/2. 

The system can be quantised in the usual way in terms 
of the bosonic annihilation and creation operators 6, 
for the probe and a, a) and for the signal cavity defined 

by[3] 




The Hamiltonian may then be written as 

H = hujpb^b + hw s a ] a + hxb^ 2 b 2 + hx a ^ 2 a 2 

+ hVx^Xb (&V 2 + fet v + 4a t a6 t 6 ) ( 6 ) 

Unlike Kumar and Divincenzo[l] we have not neglected 
the terms like b 2 a) 2 as we will choose ou p = uj s so that 
these terms are resonant. 



We now include the dissipative channels for this model 
in the usual way. The density operator for the total sys- 
tem, in the interaction picture, satisfies 

^ = -i[H T , p] + n a V[a)p + n h V[b)p (7) 

where V[L]p = LpL^ - \{pLp + ptfL) and 

Hi = h X bb ] 2 b 2 + h Xa a ] 2 a 2 + fi(e*6 + eb ] ) 

+ hVx^a OV 2 + 6 f V + 4aW6) (8) 

and K a ,Kb are the decay rates of the photon number in 
the signal and probe cavity respectively and we have in- 
cluded a resonant coherent driving of the probe cavity 
with e = sfl%Eb where \eb\ 2 is the photon flux of the 
driving field. We have also assumed that each cavity sees 
a zero temperature environment. 

In the absence of the SQUID mediated interactions the 
probe cavity will relax to a coherent state with the steady 
state amplitude 



We will chose the phase of the probe driving as a reference 
phase and set Po to be real. If we make a canonical 
transformation to the displaced picture by 

b = b + po (10) 

we can linearise the Hamiltonian, Eq. 8, in 6, b^ to obtain 

Hj = H a +4h^^p (b+P)a^a+2h^^ 2 ) 

(11) 

where the effective Hamiltonian for the signal mode alone 
is 

H a = hXa<J 2 2 + ^yJXaXbPl ^ a + h^XaXbPo {o 2 + o) 2 ) 

(12) 

which is equivalently to a parametrically driven Kerr non- 
linear cavity. This model was considered by Wielinga and 
Milburn [4]. It is equivalent to a double well system with 
a hyperbolic fixed point at the origin in phase space and 
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FIG. 1: A persistent cat Here we look at the effect of a 
including the dephasing term in addition to the the bath of 
two-photon absorbers to the ring initially in its ground en- 
ergy eigenstate. a The steady state solution for the envi- 
ronment as derived in this supplementary material and de- 
scribed by master equation Eq. (13) and (14) with T 2 = 0.2 
(or L 2 = V&2a 2 ) and T± = 0.05 (or L± = VoMa^a) and 
we have set the damping coefficient n a — 0. b We show that 
even for an environment, other than the one considered in a, 
where dephasing dominates over the the two-photon absorp- 
tion process (L 2 = a/0. 02a 2 and L± — \/0.0Sa^a) it is still 
possible for the steady state of the ring to be a cat. 



two elliptic fixed points symmetrically displaced from the 
origin. The second and third terms in Eq. (11) can be 
given a familiar interpretation. The second term is of the 
same form as the radiation pressure interaction between 
a mechanical resonator (b,b^) and a cavity field (a,af). 
The last term is equivalent to the quantum derivation 
of sub/second harmonic generation considered by Drum- 
mond et al.[5]. 

We now assume that Kb, the line width of the probe 
cavity is large, Kb » ^a^y/XaXb an d we adiabatically 
eliminate it from the dynamics. In that case from the 
point of view of the signal mode, the first term in Eq.(ll) 
looks like a fluctuating cavity detuning while the last 
arms looks like a two photon loss term. This can be ver- 
ified by explicit adiabatic elimination of the probe cavity 
field. We assume that the probe cavity, in the displaced 
picture, remains very close to its steady state of zero 
photons. The method is described in [6]. The effective 
master equation for the signal cavity is 

= -\[Ha,Ps] + T 2 V[a 2 ]p 8 + T ± V[a^a]p a + K a V[a]p s 

(13) 

where the two photon decay rate T 2 and dephasing rate 
T^ are given by 

1 2 = 

Kb 9 (14) 

1 _L = 1 2/4 = 

K b 

A peculiar feature of using SQUID coupled cavities 
is that the price paid for two-photon decay is an addi- 
tional dephasing term on the signal cavity field. Using 
the strong dependance on the steady state amplitude /?o 
in the two photon rate we can make this term domi- 
nate over the single photon decay of the signal cavity 
over the time scales of interest. In Fig. 1 we show that 
the dephasing term that is introduced in the above (un- 
damped, K a = 0) master equation has little effect on the 
Schrodinger cat nature of the steady state solution asso- 
ciated with the two-photon absorbing bath. We therefore 
believe that the discussion in the main article is in-line 
with the behaviour of realistic environments. 



[1] S. Kumar and D. Divincenzo, Phys. Rev. B 82, 014512 
(2010). 

[2] Chunqing Deng, J. M. Gambetta, and A. Lupasu, Phys 

Rev. B 82, 220505 (2010). 
[3] M. Wallquist, V.S. Shumeiko and G. Wendin, Phys. Rev. 

B 74, 224506 (2006). 
[4] B. Wielinga and G.J. Milburn, Phys. Rev. A, 48, 2494- 



2496, (1993). 

[5] P.D. Drummond, K.J. Mcneil, D.F. Walls, Optica Acta, 
28, 211 (1981). 

[6] D. H. Santamore, Hsi-Sheng Goan, G. J. Milburn, and M. 
L. Roukes Phys. Rev. A 70, 052105 (2004). 



